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Abstract 

Galilean invariant Schrodinger equations possessing nonlinear terms coupling the 
amplitude and the phase of the wave function can violate the Ehrenfest theorem. An 
example of this kind is provided. 

The example leads to the proof of the theorem: A Galilean invariant Schrodinger 
equation derived from a lagrangian density obeys the Ehrenfest theorem. The theorem 
holds for any linear or nonlinear lagrangian. 
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1 Introduction 



The linear Schrodinger equation for a pointlike object interacting with an external force 
derived from a potential U{r), 



,| = - + Uir)^ (1) 

obeys the Ehrenfest theorem[l]. The theorem is expressed as an equation for the velocity 
and an effective second law of Newton. The probability function for the averaging procedure 
is taken to be the absolute value squared of the wave function that solves the pertaining 
Schrodinger equation. 

The Ehrenfest equations read 

d . 
<v> = — <r> 
dt 



< r > 



= J d^x \ip{r,t)\'^ r 

(2) 



and 



m— < V > = F(t) 
dt ^ ^ 

F{t) = - 1 d^x \ij{r,t)\^ V Uir) (3) 

The theorem states that the quantum expectation values of position and velocity of a 
suitable quantum system obey the classical equations of motion. This observation supports 
the choice of Born[2] of the amplitude squared as a probability weight factor. According to 
the present interpretation of quantum mechanics [3], the outcome of measurements of the 
motion of identically prepared pointlike systems, their location and speed when subjected to 
so-called 'external' forces, will be, on average, the same as the sharp values of their classical 
counterparts. 

In classical deterministic mechanics, the sources of randomness are twofold, statistical 
fluctuations and extreme sensitivity to initial conditions. Statistical fluctuations can be 
minimized by bringing the system in contact with a reservoir at very low temperature, 
near absolute zero. The sensitivity to initial conditions causing chaotic behavior are not 
omnipresent. They arise for nonlinear systems. Hence, the outcome of the measurement of 
a point particle is in principle, predetermined. Eqs.(l) tell us that besides the statistical 
fluctuations, quantum systems possess an extra source of indeterminacy, that is regulated 
in a very definite way by the complex wave function. This is in some sense better for the 
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experimenter than classical thermal fluctuations that are hidden completely. Even in the 
limit of zero temperature, at which the classical fluctuations vanish and a deterministic 
outcome is expected, the macroscopic variables position and velocity are still fuzzy. The 
fuzziness is prescribed by the Schrodinger equation. Hence, only averages obey deterministic 
equations. 

Ehrenfest theorem can be extended to many point particle systems without difficulty. 
The theorem embodies Bohr's correspondence principle. As such, the dynamics dictated by 
eq.(l) ought to be independent of Planck's quantum of action h. A failure of the ehrenfest 
theorem would carry on to the behavior of quantum particles in external fields such as 
gravitational fields. A system that is assumed to obey a nonlinear Schrodinger equation of 
the kind described below, displays an inherent nonclassical behavior, even on the average. 

In section 2 we single out the kind of nonlinearities that violate with Ehrenfest theorem. 
In section 3 we prove a theorem that connects Galilean invariance, and the existence of 
a lagrangian density whose Euler- Lagrange equation is the Schrodinger equation, to the 
fulfillment of the Ehrenfest theorem. The proof hinges upon the existence of a lagrangian 
density. The existence of this density guarantees the existence of a conserved energy and 
a conserved linear momentum for the free particle, even when the particle is allowed to 
self- interact. Galilean invariance then inforces a connection between the Ehrenfest average 
velocity and the conserved momentum. 
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2 Ehrenfest theorem breaking 



The linear Schrodinger equation for a pointlike quantum particle whose configuration wave 
function is in interaction with an external potential U{r), can be derived from the 

lagrangian 

The lagrangian is a real scalar. A global phase transformation on the wave function 
leaves the lagrangian invariant. The Noether current for the symmetry becomes [4] 

jo ^ ip* i; 

J = (vr V' - r) (5) 

2 m V / 

Born [2] interpreted this current as a probability current. Jq, the probability density, and, 
j, the probability flux. For waves whose amplitude decreases fast enough at infinity, current 
conservation 



f.v.J^o (a) 

function to unit norm, the modified norm is in Born's view the total probability to have the 
particle anywhere in space. 

Ehrenfest[l], before Born, found that there is a way to define a privileged location whose 
dynamics follows the classical equations of motion. After Born's interpretation, it becomes 
clear that the definition 

< r >^ J d^x r (7) 

is a sound one for the averaged position variable or packet centroid. Derivation with 
respect to time of eq.(7), application of eq.(6), and integration by parts for well behaved 
wave functions yields eq.(2) with 

< v(t) > = j d^x^ (8) 

Further derivation of the latter and use of the Schrodinger equation of eq.(l), yields the 
eq.(3). 
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Consider now nonlinear Schrodinger equations. Such equations as the Gross- Pitaevskii 
equation[5] are frequently used for the description of Bose-Einstein condensates (BEC) of 
Alkali gases. The equation is an effective equation, derivable from a field theory by taking 
condensate expectation values, and it has found quite success in the description of the BEC 
setup. Starting from a linear field theory of particles, it is possible to arrive at a nonlinear 
Schrodinger equation in which the interactions arc of zero range to lowest order. In this 
sense the effective theory is merely an approximation to the actual complicated dynamics 
ruling the behavior of the particles in interaction. As an effective theory, it has found a fair 
amount of predictive power. It is logical then, to ask the question whether this effective 
theory still obeys the Ehrenfest theorem. If the answer is in the affirmative, then, the 
averaging procedure that lead to the Schrodinger equation for the condensate did not break 
the connection to classical physics. If not, then, either this breaking is a falsification due to 
the approximations, or else, it really refiects an observable feature of these systems. In our 
opinion, only experiment can answer this question. 

Another widespread application of nonlinear Schrodinger equations is that of density 
functional phenomenological models. In these models, the self-interaction of the particles, 
is nonlinear in the density. This approach yields statical and dynamical properties of sys- 
tems, ranging from electron transport in solids[7], electronic excitations [8], soft condensed 
matter[9] phase transitions in liquid crystals[10], phonons in solids [11, 12], colloids[13], liq- 
uids and nuclei [14], atoms and molecules [15], quantum dots[16], etc. In these theories, the 
Schrodinger equation becomes 



^^ = - + [0{m + U{rM (9) 

where Od-i/"!^) is a nonlinear and sometimes nonlocal[6] functional of the density p = \^Jj\'^. 
Calculating the time derivative of the velocity of eq.(8), with the nonlinear term of eq.(9) 
included, we find 



m — < v(t) > = - fd^x p(r)- V [U{r) + 0{p)] (10) 
Using the trivial property 



I d^'x p V[0{p)\ = (11) 

we arrive at the same Ehrenfest equation as that of eq.(2, 3). Nonlinear Schrodinger 
equations with nonlinearities such as that of eq.(9) obey the Ehrenfest relations. Perhaps 
this is expected as the nonlinear terms are essentially classical, depending on the density p 
and not on the wave function itself. If the external interaction is itself nonlinear such as 
U{p,x), then the Ehrenfest therem is still satisfied, but this time with a nonlinear force ^ 

^This is analogous to the results found in soliton models. [17] 
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f{t) = -jd^x\^{r)\'^VU{p,r) (12) 

In the following we will show that nonlinear Schrodingcr equations of a different type, 
for which the nonlinearity depends on the phase also, do not obey the Ehrenfest equations. 
The trivial property of eq.(ll) is the clue to the construction of the desired equations. The 
vanishing of eq.(ll) is due to the dependence of 0(p) solely on p. It is easy to show, that an 
operator that depends on gradients of the phase for instance, will yield extra terms to the 
Ehrenfest equations. 

Consider the Schrodinger equation proposed by Doebner and Goldin[18], without the 
diffusive term ^ 



(13) 

where A is a coupling constant and j is defined in eq.(5). Current conservation still holds 
as in eq.(6). The definition of the velocity is still given by eq.(2). Equation (13 is Galilean 
invariant. Galilean invariance is effected by means of the transformations [19, 4, 20] 



r r — S-v t 

(j) ^ {5wf - mSv ■ r (14) 

with 6v a constant velocity parameter. Under these substitutions, the free Schrodinger 
equation is invariant , while 



-4 + (15) 



Eq.(13) is then Galilean invariant. However the Ehrenfest theorem is not satisfied 



m ^ < v(t) >= m + X J d'x V|^|2 • v(^) (16) 

The identity (11) is of no use now. The last term in eq.(16) depends on derivatives of 
the amplitude and derivatives of the phase of the wave function, and not on p solely as in 

eq.(ll). 

sizeable body of literature on related topics may be traced by citations to the ref.[18] 
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The failure to satisfy the Ehrenfest theorem may by traced back to the absence of a 
lagrangian density for which eq.(13) is its Euler-Lagrange equation of motion. This remark 
will be put in the form of a theorem in the next section. 
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3 Ehrenfest theorem as a consequence of Galilean in- 
variance 



The proof of the theorem proceeds by the following steps. The Schrodinger lagrangian for a 
free particle including self-interactions of any nonlinear nature, but no explicit dependence 
on the space or time coordinates is introduced. The corresponding action is then invariant 
under spatial coordinate translations. By Noether's theorem there arises a conserved current 
and the physical law of conservation of linear momentum. The lagrangian is also demanded 
to be a real scalar. It depends on the phase of the wave function only through its derivatives. 
Phase transformations will then induce the law of conservation of probability identified as 
the modulus squared of the wave function. 

Galilean invariance of the lagrangian then determines a connection between the proba- 
bility current and the linear momentum. Finally, this connection insures the validity of the 
Ehrenfest theorem. 

Consider a lagrangian density 



C{ilj,ipt,'4'i) = C{St,Si,R,Rt,Ri) 

= Re^^ (17) 

where suffixes denote partial derivatives, with respect to time t and with respect to the 
coordinate Xi. As we are in the nonrelativistic domain, we do not differentiate between 
covariant and contravariant indices. 

Space translations are generated by the transformations f — > r -|- e, with e a constant 
parameter. As the lagrangian is independent of the coordinates, it is invariant under this 
transformation. When the parameter e is promoted to be space-time dependent [4] , the law 
of conservation of linear momentum appears as a condition on the invariance of the action 
for arbitrary e. 

The variation of the action reads 



5A 



d d X dt C 



j d^x dt 



CV.e--e,.VS-^e,.VR 
^^-VR.e, 



dRi 



(18) 



The invariance of the action then engenders the law of conservation of linear momentum 



= 



dt 



dxi 



(19) 
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where 



p = 


dC 
dSt 


VS- 


oRt 




dC 

d ipt 


Vip- 






SijC - 


dC 


B 

dRi 




SijC - 


dC 





As the lagrangian (17) is independent of S, a transformation of the form 5 — > + t), 
will generate the law of conservation of probability. Straightforward manipulations lead to 



(21) 

where 



Jo 



— * 

J 



In both equations (20,22) we do not specify the dynamics, by replacing the lagrangian 
with the free Schrodinger lagrangian. The only conditions imposed are, the very existence 
of a lagrangian that is a real scalar dependent on a complex wave function and indepen- 
dent of spatial coordinates except through the wave function and its derivatives. We also 
assumed that there do not appear derivatives of higher order, although this is not an es- 
sential ingredient. Both eqs. (19,21) may be obtained also by means of the Euler-Lagrange 
equations [21, 22]. It is easy to check that they are mere identities when the Euler-Lagrange 
equations are used. 

The connection to Ehrenfest theorem now proceeds through the demand of Galilean 
invariance of the Schrodinger equation or covariance of the action for finite Galilean boosts. 



"dSt 
dC 



- I 



dC 



dVS 

9j0 . dC . 
I — — ip — I — — ip 



(22) 
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Covariance, and not invariance, as a boost modifies the kinetic energy, and consequently 
changes the action. Applying an infinitesimal galilean transformation to the lagrangian L 
as specified in eqs.(14), i.e. dropping the term quadratic in the velocity, the variation of the 
lagrangian becomes ^ 

dSt dRt dVS 

— 5v • p — mSv • J (23) 

where we have used the definitions of p and J in eqs. (20,22). For the action to be 
covariant, which is equivalent to being invariant for infinitesimal Galilean boosts, eq.(23) 
implies 



P = mJ + V/(f,i) (24) 

Galilean invariance requires the probability fiux to differ from the conserved linear mo- 
mentum by at most a total divergence. For all practical purposes both vectors are one and 
the same up to the factor of the mass. 

The existence of a real scalar lagrangian, even including nonlinear self-interaction terms, 
and the requirement of independence on coordinates for the free self-interacting is neces- 
sary and sufficient for the equivalence between m J and p. Necessary, as we can proceed 
backwards from eq. (24) and reconstruct Galilean invariance. 

We are now ready to show the validity of Ehrenfest theorem for a generic lagrangian as 
in cq.(17). 

We add a scalar potential term = —U{r) \^\'^ to the lagrangian. Use of different types 
of potentials, such as velocity dependent vector potentials, do not change the conclusions. 
Moreover, nonhnear potentials will generate forces as in eq.(12). 

The lagrangian is no longer invariant under translations, but the action still is. 

Equation (19) is now modified to 

Integrating over space for asymptotically vanishing wave functions, we find 



fj d'x p,iv,t) = - I d'x^^\ij\' (26) 



This is the second law of Newton for the 'field' momentum p. It bears no relation to the 
Ehrenfest theorem at this stage. 

^This method is straightforward for our purposes, use of the Noether method yields a tautology. 
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Consider now the definition of < f > of eq.(2). Using eq.(21) we find 




d 



J dt 



-J d^xvV-i 
j d^xJ 



(27) 



Using eqs. (24,26) we find the second law of Newton for the centroid of the packet < r >, 
the Ehrenfest theorem 



We have shown that the demand of Gahlean invariance on a scalar real local lagrangian 
pertaining to a Schrodinger field, linear or nonlinear in the wave functions, implies the 

validity of the Ehrenfest theorem. The same does not hold for Schrodinger equations that 
are not derivable from an action principle. Among the terms that fulfill the requirements of 
Galilean invariance, and can serve as viable nonlinear self inter-actions, we can mention the 

— * — » 

function g{q), with q — Rt + j • 'Vln{R). 




= < F > 



(28) 



4 Conclusions 
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